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FINAL REPORT 


Work on this project began with a laboratory investigation of ion reaction schemes 
that would allow the near simultaneous measurement of NH3, DMSO, and DMS 0 2 with 
a single initial reactant ion. If the proper reactant ion could be found, it would make 
possible fast and sensitive measurements of the above three compounds in remote marine 
environments, which could add much to our understanding of the DMS oxidation process 
and the formation of new aerosol particles. 

Several different reactant ions have been tested as a possible means for allowing 
selected ion chemical ionization mass spectrometric (SICIMS) measurements of NH3, 
DMSO and DMS 0 2 . All of the tested ions have a sufficiently low proton affinity that 
they can directly proton exchange with the above 3 compounds of interest. In all but one 
case, however, their proton affinity is still sufficiently high to minimize reactions with 
most other atmospheric constituents. Acetone was used to produce C3H70 + ions which 

will proton exchange with NH 3 to form NH4 for the measurement of ammonia. 

However, acetone cluster ions C 3 H 7 0 + (C 3 H 6 0 ) also formed rapidly, and do not appear 
to react quickly with NH3 thus reducing the attractiveness of using these ions for the 
above measurements. Protonated DMS was thought to be an attractive reactant ion 
because DMS can reach concentration of hundreds of ppt in areas where the proposed 
measurement techniques might be used. If CH 3 SCH 3 H + ions were used, enhanced DMS 
concentrations should have little to no effect on the ion spectrum because the 
concentration of DMS used to generate the initial DMS reactant ions would be in excess 



of that present in the atmosphere. It is difficult to find very pure DMS, however, and 
compounds such as DMSO are common contaminants. Butanol ions were also examined 
as a possible initial reactant ion for the proposed measurement scheme. There appeared 
to be some potential contaminants associated with this compound also, but they may not 
be coming from the original butanol source itself. Thus, this compound needs to be re- 
examined. The H 3 0 + ion was also examined as a means for measuring NH 3 
concentrations. This appeared quite promising for remote marine environments, but a 
higher proton affinity compound with less potential interferences was still desired. 

Protonated ethanol ions were also tested and appear to offer promising new 
measurement capabilities for a number of compounds. The proton affinity of ethanol is 
fairly low (188 kcal/mol), 1 so it should react readily with compounds having proton 
affinities near or above 200 kcal/mol. While it clusters with itself like acetone, the clusters 
still appear to proton exchange quickly with compounds such as NH 3 . Protonated 
ethanol ions were used in a preliminary continental field study which demonstrated the 
usefulness of this ion in conjunction with the SICIMS technique for the measurement of 
gas phase NH 3 . This field study will be discussed later in this report. 

Much progress has also been made on our nucleation studies. This effort began 
with the design and construction of a new type of selected ion chemical ionization source: 
one in which relatively high concentrations of a specific ions could be generated and 
subsequently injected into a separate drift/reaction region. The drift/reaction region is also 
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unique in that it allows reaction times to be varied from less than 1 ms to more that 1 00 
ms in an atmospheric pressure gas flow. 

This wide range of reaction times is important because it provides a means of 
separating ion induced molecular cluster formation from ion clusters which result from 
proton exchange with a pre-existing neutral H 2 S0 4 dimer or trimer. In the absence of 

neutral H 2 S0 4 dimers or trimers, HS0 4 (H 2 S0 4 ) n ions (with n > 0) will still form in 
the presence of gas phase H 2 S0 4 by successively incorporating H 2 S0 4 molecules onto 

HS0 4 ions until some steady state value of n is reached. The latter HS0 4 (H 2 S0 4 ) n 
ions (which will be called ion-induced clusters) are indistinguishable from those resulting 
from H 9 S 0 4 dimers or trimers which loose a proton to some initial reactant ion. 1 hus, 

HS0 4 • (H 2 S0 4 ) n ions which are being used to investigate neutral H 2 S0 4 cluster 
formation have at least two sources: 


NOj + H 2 S0 4 -» HS0 4 + HNO 3 (la) 

HS0 4 +H 2 S0 4 + nw HS0 4 H 2 S0 4 +m (lb) 

HS0 4 (H 2 SO) n l +H 2 S0 4 + m-» HS0 4 -(H 2 S0 4 ) n +m (lc-n) 

and NO 3 + (H 2 S0 4 ) n — > HS0 4 (HS0 4 ) n _, + HNO 3 (2) 
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Reaction 1 is a multi-step reaction in which the ratio of jHSC >4 • (H 2 S0 4 ) n J/JhSO^] is 
dependent on [H 2 S0 4 ] and reaction time. In contrast, reaction 2 is simply a proton 
exchange reaction between any of the (HS0 4 ) n clusters and NO 3 , and the ratio of 
[HS0 4 (H 2 S0 4 ) n j]/[hS 0 4 ] should not decrease rapidly as reaction time is shortened. 
Thus, if the HS0 4 • H 2 S0 4 ]/|hS0 4 ratio is plotted as a function of time for short 

reaction times as shown in Figure 1, a linear plot would be expected for the latter ratio if 
HS0 4 • H 2 S0 4 were formed solely by an ion induced process. Figure 1 shows data 
obtained using the recently constructed 1-100 ms ion source/reactor apparatus. The 
approximately linear time dependence is expected at the relatively long reaction times 
shown, because under such conditions ion-induced cluster formation is known to be 
important. At even long times, the ratio of [hS0 4 • H 2 S0 4 ]/[HS0 4 ] must eventually 
become time independent at some equilibrium concentration ratio. As the reaction time is 
decreased, however, if (H 2 S0 4 ) 2 exists and reacts with NOJ (which it is expected to), 
then the slope of the line in Figure 1 is expected to approach a constant in which the ratio 
of [(H 2 S0 4 ) 2 ]/[H 2 S0 4 ] is represented by |hS0 4 • H 2 S0 4 J/[hS0 4 times some 

relative reaction rate constant for N 03 /(H 2 S 0 4 ) 2 and NC> 3 /H 2 S 0 4 . Therefore much of 
our efforts have gone into trying to measure this ratio under near ambient conditions and 
at the shortest possible reaction times. The near linear dependence of the ratio (even 
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H 2 S0 4 monomer / dimer signal 



down to the 1 ms time scale) suggests that less than 1 % of the gas phase H 2 SO 4 
molecules are dimers, even at the extreme cases of -20°C and [H 2 SO 4 ] — 100 pptv shown 
in Figure 1 . Thus, at temperature and [H 2 S0 4 ] more typical of previous field studies 
where nucleation events have been observed, the percentage of (H 2 S0 4 ) 2 clusters 

present would probably be expected to be orders of magnitude lower still. While this 
does not rule out the role of H 2 S0 4 dimers and trimers in particle nucleation it does make 
it seem far less likely that a sufficient number of these clusters exist, and that they will 
react fast enough with NH 3 to control tropospheric nucleation. As a consequence of 
these laboratory studies, the most likely tropospheric nucleation mechanism would now 
appear to be the gas phase reaction of NH 3 with H 2 S0 4 (rather than with H 2 S0 4 dimers) 

to form ammonium bisulfate molecules and their subsequent formation of clusters . 2 ' 4 The 
present study has caused us to reconsider the relevance of a previous experimental 
laboratory investigation aimed at determining the potential interference of NH 3 in H 2 S0 4 
measurements. This previous data suggested that while NH 3 does not react quickly with 
H 2 S0 4 in the gas phase, it may, in fact, react slowly. Determining the extent to which 
NH 3 and H 2 S0 4 react to form ammonium bisulfate has, thus, become a central focus of 
the present study. 

A flow tube and injector system has recently been interfaced with the new ion 
source and SICIMS instrument in an attempt to measure the rate coefficient for the gas 
phase H 2 S0 4 /NH 3 reaction. This is an extremely difficult measurement, for several 



reasons both H 2 S0 4 and NH 3 can be rapidly deposited onto the flow tube walls, they can 
also react with each other on the walls, and the gas phase reaction of interest is expected 
to be very slow. These measurements are presently under way. We hope to complete 
them over the next few months and to subsequently look for the molecular cluster 
products expected from this reaction. This project is divided into two parts because of an 
administrative fund change at NASA. (While this is a final report for the first portion of 
this project the overall research project itself is expected to continue for an additional year 
as originally proposed. It is during this final year that the latter reaction would be 
studied.) 

Other areas of progress include measurements of the accommodation coefficient of 
H 2 S0 4 on ammonium sulfate, sodium chloride and other particles surfaces. These 
measurements were made in our laboratory in conjunction with colleagues from the 
University of Minnesota. The accommodation coefficient for H 2 S0 4 is important not 
only for understanding the growth of CN into CCN, but it is also used to calculate the 
flux of nucleating particles from measured ultrafme particle distributions. Since literature 
values of this coefficient vary by nearly 2 orders of magnitude (0.02-1 .0), and we had the 
ability to make the first direct measurement of these parameters, it seemed like a 
worthwhile project. The accommodation coefficients were found to be 0.73 ± 0.21 and 
0.80 ± 0.23 for ammonium sulfate and sodium chloride respectively. The details of this 
very successful measurement have recently been published and a copy of the paper is 
included as Appendix A. 
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Progress has also been made on testing new ion reaction schemes for measuring 

NH 3 in the field. In the first year of this grant, several ions were tested in the laboratory 

* 

as possible candidates for the reactant ion to be used for measuring NH 3 in the 
atmosphere. This year, an approximately 1 month long local field study was conducted, 
in which two reactant ions were tested for use in measurements of ambient NH 3 and 
various hydrocarbons by SICIMS. One of the compounds tested was acetone, which was 
previously used in some of the past year’s laboratory tests. A second compound used to 
generate reactant ions was ethanol, and the latter results look very encouraging. Both 
compounds readily form ions by attaching a proton, and both also readily form clusters 
with themselves. In the case of ethanol the most common cluster ion C 2 H 70 + C 2 HgO 

still appears to react quickly with NH 3 via proton exchange to form NH 4 and 
NH 4 • C 2 H 6 0. Figure 2 shows the initial ethanol reactant ion distribution in ambient air 
which is dominated by C 2 H 7 0 + • C 2 H 6 0, but which also contains a small peak of the 
core ion C 2 H 7 0 + . The NH 4 and NH 4 • C 2 H 6 0 peaks are a result of reactions of 
C 2 H 7 0 + ■ (C 2 H 6 0) n for n=0,l with ambient NH 3 at a field site a few km north of the 

NCAR Mesa Laboratory Building at the edge of the foothills of the Rocky Mountains. 
Air flow at this site is typically from the west and is relatively clean during the morning 
hours. At some point during late morning or early afternoon, the flow will typically 
switch to a more easterly (or sometimes northerly and southerly) component, which 
brings up air from the city of Boulder, Denver, or farming areas to the east. Figures 2 and 
3 show a measured ion spectrum and an initial test of the new SICIMS ammonia 
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Typical ion spectrum measured during a recent field study showing 
nearly exclusively reactant and production peaks. 
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measurement technique used at this site. Figure 3 also shows the wind direction during 
this same time period. As is typical of this site, prevailing westerlies bring down 
relatively clean air from the mountains which, as seen from Figure 3, contains little NI I 3 . 
In early afternoon at about 14:00 hours, the wind switches direction abruptly and its 
velocity drops. The air mass present after 14:00 is probably from urban and farm areas 
and would be expected to contain far more NH 3 , which agrees well with the observations 
shown in Figure 3. Thus far, only very preliminary tests of this new technique have been 
performed, but the results look quite promising, particularly since the area in which these 
measurements were made contains large amounts of hydrocarbons which could have had a 
much larger than observed interference. For measurements in marine environments, 
C 2 H 70 + can probably be used to measure NH 3 , and DMSO with the same ion source. 

The present development and testing of new SICIMS measurement schemes is 
also expected to continue throughout the next year. While ions such as C 2 H 7 0 + show 
much promise for making measurements of compounds such as NH 3 , and DMSO, it is 
quite possible that DMS and DMS0 2 can also be measured with the same ion source. 
Thus, these and similar reaction schemes for other DMS oxidation products will continue 
to be a major focus of the second phase of this study. 
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APPENDIX A 
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Measurements of the H2SO4 mass accommodation coefficient 
onto polydisperse aerosol 

A. Jefferson , 1 F. L. Eisele , 1 ’ 2 P. J. Ziemann , 3 R. J. Weber , 4 J. J. Marti , 5 
and P. H. McMurry 6 

Abstract. The loss rate of H2SO4 vapor onto submicron particles was measured for three 
different particle substrates. The experimental technique involved direct flow tube measurements 
of H 2 SO 4 decay rates onto a polydisperse aerosol using chemical ionization mass spectroscopic 
detection. The aerosols of this study were partially hydrated crystalline salts with diameters in 
the size range of 20 to 400 nm. The mass accommodation coefficients, a, were calculated from 
the first-order rate constants for H 2 SO 4 loss to be 0.73 £ 0.21 and 0.79 ± 0.23 for loss onto 
(NH 4 ) 2 S 04 and NaCl, respectively. Measurements of the loss rate of H 2 SO 4 onto a NaCl 
aerosol coated with stearic acid resulted in lower mass accommodation coefficients with values of 
0.31 and 0.19 for aerosol with high and low stearic acid coverage, respectively. The observed 
decrease in a on an aerosol with a hydrocarbon coating suggests that aerosol composition is a 
key factor in H 2 SO 4 adsorption on to a particle surface. 


Introduction 

Aerosols have a significant impact on climate forcing by 
their ability to scatter solar radiation and increase cloud 
albedo, which result in a cooling of the lower troposphere 
[Permer et al. , 1997; Wigley and Paper, 1992], The magnitude 
of this effect depends in part on the aerosol microphysical 
properties, particularly the particle composition, size, and 
number density [Schwartz, 1996]. Sulfates, along with nitrates 
and organics, form the primary constituents of submicron 
condensation nuclei (CN), which scatter shortwave radiation 
and can grow to form cloud droplets. The formation and initial 
growth of CN are thought to proceed through bimolecular 
nucleation ofH 2 0 and H 2 S0 4 and subsequent uptake of gas 
phase H 2 S0 4 onto the particle surface [Kriedenweis et al ., 
1991], An unknown parameter needed to quantify CN 
nucleation and growth rates is the mass accommodation 
coefficient of H 2 S0 4 onto particles. 

Current models of aerosol nucleation rates are especially 
sensitive to the H 2 S0 4 mass accommodation coefficient, a. 
Estimates of a range between 0.02 and 1.0. In their mode! of 
bimolecular nucleation rates, Pandis et al. [1994] calculated a 
45% reduction in the H 2 S0 4 -H 2 0 nucleation rate when a was 
changed from 0.02 to 0.05. Nucleation rates modeled by 
Kulmala et al. [1995] show a similar strong dependence on a. 
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In laboratory measurements of aqueous H 2 S0 4 particle 
growth. Van Dingenen and Raes [1991] determined a value 
for a between 0.02 and 0.09. Alternatively, two field studies, 
which measured steady state H 2 S0 4 production rates [Eisele 
and Tanner , 1993; Weber et al ., 1997], found that values for a 
of 0.5 and 1.0, respectively, gave good agreement with their 
data. 

Submicron aerosol composition may vary' between remote 
or polluted and marine or continental environments. Aerosols 
from remote marine environments were found to contain mostly 
S0 4 2 ', CH3SO3-, and NH 4 + ions, while larger sea salt aerosols 
also contain Na + and CT ions [Huebert et al. , 1996, Savoie et 
al ., 1993]. In addition, aerosols found in continental and 
anthropogenically influenced air masses may contain 
considerable amounts of organics [Saxena and Hildemann , 
1996; Novakov and Penner , 1993]. As the substrate 

composition varies, so too may the uptake of H 2 S0 4 . As an 
acid, H 2 S0 4 would seem likely to bond most strongly with 
substrates containing polar or basic binding sites. 

Here we present measurements of the mass accommodation 
coefficient of H 2 S0 4 onto aerosols composed of (NH 4 ) 2 S0 4 , 
NaCl, and NaCl coated with stearic acid. This study is the first 
measurement of the gas phase H 2 S0 4 uptake rate onto aerosols 
that directly measured changes in the gas phase H 2 S0 4 
concentration. The aerosol uptake experiments in this study 
use the sensitive technique of selected ion chemical ionization 
mass spectrometry (SI/CIMS) to measure the removal rate of 
H 2 S0 4 by aerosols. In this study, H 2 S0 4 concentrations in 
the range of 0.1-10 parts per trillion by volume were measured 
at a relative humidity at or below 11%. Measurements at a 
higher relative humidity, representative of the lower 
troposphere, were hindered by the low equilibrium vapor 
pressure of H 2 S0 4 above the acid solution source used in 
these experiments. Possible changes in the H 2 S0 4 mass 
accommodation coefficient with ambient conditions such as 
relative humidity and aerosol composition are discussed. 

Experimental Method 

Polydisperse aerosols were generated by atomizing 
aqueous solutions of either (NH 4 ) 2 S0 4 or NaCl. Filtered 
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nitrogen was used as a earner gas for the aerosol flow at a flow 
rate of about 1.24 L min" 1 . The aerosol droplets passed 
through a silica gel diffusion dryer before joining the 
humidified nitrogen earner flow at the entrance of the flow 
tube injector. Particle sizes ranged in diameter from 20 to 400 
nm. In the experiments involving stearic acid coated aerosols, 
NaCl particles were coated with stearic acid by passing the 
aerosol flow through a heated 6 inch long glass tube 
containing both solid stearic acid and its vapor. The thickness 
of the acid coating on the particles was controlled by 
adjusting the temperature of the glass tube and thus the stearic 
acid vapor pressure. The tube was heated to 104 and 88 C for 
high and low particle coatings of stearic acid, respectively. 
The aerosol flow made up about 5% of the total gas flow. 

Aerosol surface areas were determined from particle size 
distributions measured by using a scanning electrical 
mobility spectrometer (SEMS) [Wang and Flagan, 1 990], The 
instrument consists of a particle charger, a differential mobility 
analyzer (DMA), a condensation nucleus counter (CNC), and 
SEMS software for rapid scanning of the panicle size 
distribution. Sampled particles were charged to near- 
Boltzmann equilibrium by a 210 Po source. The particles then 
passed through a DMA, which separated particles by their 
electrical mobility, allowing only those with mobilities in a 
narrow selected range to exit to the CNC detector. The particle 
mobility distribution was measured by scanning through the 
particle size range, typically 20 nm to 400 nm in diameter. For 
spherical particles, whose size is precisely defined by the 
diameter, the particle size distribution can be accurately 
calculated from the mobility distribution. 

In our experiments, aerosol was sampled from the 
downstream end of the flow tube reactor and diluted with 
filtered, dry nitrogen (-1:6) in order to maintain particle 
concentrations in the instrument range. The aerosol flow rate 
into and out of the DMA was 1.0 L min' 1 , the DMA sheath air 
flow rate was 7.3 L min"^ of filtered room air, and the CNC flow 
rate was 1.5 L min' 1 , of which 1.0 L min' 1 was aerosol from the 
DMA, the remainder being filtered room air. The use of filtered 
room air for the DMA sheath flow should have no effect on the 
measured particle size distribution, since the ambient relative 
humidity was well below the deliquescence point of both 
NaCl and (NH 4 ) 2 S0 4 . 

The particle size distributions used for our calculations of 
aerosol surface areas were corrected for particle charging and 
losses in the DMA by using the method of Hagan and Alofs 
[1983]. No correction for CNC counting efficiency was 
required, since the aerosol surface area was associated with 
sizes that are detected with near 100% efficiency. The flow was 
laminar in the flow tube and the tube leading to the DMA, 
making particle losses by diffusion negligible. Aerosol surface 
areas were calculated by assuming that the particles were 
spherical, which is approximately true for (NH 4 ) 2 S0 4 [Kelly 
and Me Murry, 1992] although not for NaCl. Nonetheless, 
because NaCl particles are apparently cubic [Lodge and Tufts , 
1955] and the ratio of the side length of a cube to its mobility 
diameter is 0.746 [Dahneke, 1973], the surface area for a sphere 
and a cube with the same mobility diameters differ only by 
-6%. Typical aerosol surface areas per cubic centimeter for 
these experiments were 600 }im : cm* 3 . 

Flow Tube Apparatus 

The reaction tube for measurement of kinetic decay rates 
consisted of a jacketed Pyrcx flow tube, 80 cm in length with a 


7.6 cm inside diameter. A movable, 2.5 cm diameter Pyrex 
injector was mounted inside the flow tube. The injector 
diameter reduces to a small orifice about 5 cm before its exit. 
After this constriction the injector flares to nearly the full 
diameter of the flow tube and has a stainless steel, fine mesh 
screen attached on the end. The H 2 S0 4 vapor source consisted 
of a small cap containing a concentrated H 2 S0 4 solution 
which rested on top of the injector end screen. The injector 
design ensured that the aerosol flow turbulently mixed with 
H 9 S0 4 vapor in this region between the constriction and 
injector exit and that the total gas flow was evenly distributed 
across the flow tube diameter. A schematic diagram of the 
experimental setup is given in Figure 1. 

The H 2 S0 4 concentration in the flow gas depended on the 
equilibrium vapor pressure of H 2 S0 4 above the acid solution 
source and thus varied strongly with the relative humidity of 
the gas flow. At relative humidities above 10% the H 2 S0 4 
vapor pressures of sulfuric acid solutions drop below 10* 9 torr 
[Ayers et al, 1980; Marti et al. , 1997], which resulted in 
H 2 S0 4 concentrations too low to measure an accurate decay. 
The present measurements were limited to a relative humidity 
between 1 1% and 0.3%. 

The temperatures of both the flow tube and the humidified 
carrier gas were maintained at a constant 25°C. A humidified 
airstream was generated by flowing N 2 gas through 
microporous polytetrafluoroethylene (PTFE) tubing immersed 
in deionized and distilled water. This airstream then flowed 
through a heat exchanger before combining with the aerosol 
flow in the injector tube entrance. Water from a temperature- 
controlled bath was circulated through both the external flow 
tube jacket and the gas heat exchanger. The gas temperature 
was measured by a thermocouple at the flow tube exit. The 
relative humidity of the carrier gas was measured at the flow 
tube exit by a dew point monitor. 



QNS 


Figure 1 Diagram of flow' tube kinetic reaction cell, aerosol 
generation and detection, and CIMS instrument. 
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CIMS Detection of H 2 S0 4 

The decay in the H 2 S0 4 concentration as a function of 
injector position was measured by SI/CIMS. This technique 
involved sampling the gas from the center of the flow tube into 
a chemical ionization source. After merging with the ions from 
this source, the gas flowed through a reaction region, where 
H 2 S0 4 in the flow gas underwent a proton exchange reaction 
with N0 3 ' and was converted to HS0 4 *. The reactant and 
product ions were then focused into a mass spectrometer and 
detected. The relative H 2 S0 4 concentration is given as the 
ratio of the HS0 4 7N0 3 - ion signals. Two steps were taken to 
assure that hydrates of H 2 S0 4 were detected along with the 
monomer. Before entering the mass spectrometer the 
ion-containing gas was directed through a flow of dry 
nitrogen gas, which acted to dehydrate the ions. Next the ions 
passed through a small aperture into a collision dissociation 
region. A voltage held across the vacuum chamber accelerated 
the ions and induced collisions with the background nitrogen 
gas strong enough to dissociate any remaining ion clusters. 
Thus hydrates of HS0 4 ‘ were converted to the monomer sulfate 
ion- and detected with the same efficiency. Further details of 
the CIMS instrument can be found in prior publications 
[Eisele and Tanner , 1991, 1993; Tanner and Eisele , 1 995; 
Tanner et al . , 1996], 


tube, of time (in seconds) = 0.06 16*. The difference between 
this flow rate and the measured value of 0.0639* given in 
equation (1) is due to difftisional broadening of the time 
profile. The time offset of 0.3715 s in equation (1) reflects the 
time needed for the ions to travel from a relative zero position 
on the flow tube to registering an ion signal on the computer. 
The small difference in the ideal laminar flow rate and the 
measured ionic flow rate indicates that the difference between 
the radial and axial diffusion of the ions and neutral molecules 
is small for the low ion concentrations of these experiments 
and has a negligible effect on the measured flow times. The 
exponential fits of the experimental decays used a value for the 
reaction time given by equation (1). 


Reaction Kinetics and Mass Transport 

The loss of H 2 S0 4 onto an aerosol surface follows pseudo 
first-order kinetics, provided the aerosol number density and 
uptake ability remain constant during the course of a single 
measurement. Changes in the H 2 S0 4 concentration with 
injector position can be expressed as in the following 
equation; 

4H 2 S0 4 ] = _^H 2S04] + Jfcre[ H 2 S0 4 ] s (2) 


Flow Time Measurements 

At a flow rate of 23.5 L min' 1 and a Reynolds number of 
650 typical of the present experiment, the gas in the flow tube 
shown in Figure 1 is expected to develop a laminar flow 
profile with some radial and axial diffusion. However, to 
promote a short ion-molecule reaction time, the flow in the 
CIMS sample inlet was considerably faster than that in the 
main flow tube, and as a result the flow accelerated, permitting 
a radial convergence of the gas in the flow tube into the sample 
inlet. Only a central portion of the full flow tube was sampled. 
Rather than assume laminar flow conditions, we measured the 
flow time directly, using a pulsed ion source. A corona 
discharge ion source was constructed by running a wire down 
the center of the injector to a distance within 2-3 cm of the end 
metal screen, which was held at electrical ground so that the 
flow tube remained a near field free region. A counter card 
triggered a 5 kV, 10 ms wide electrical pulse to induce the 
corona discharge and also initiate computer data collection. 
Data was collected in 10 s cycles with 0.2 s time steps. The 
flow times were calculated by using a weighted average of the 
signal counts at each time step are best fit by the equation 

time (in seconds) = 0.0639* + 0.3715 (1) 

where * is the injector position in centimeters from the top of 
the CIMS sample inlet tube. 

Unlike a neutral gas, a positively charged ion flow is 
subject to space charge repulsions, which accelerate the ions, 
thereby increasing the diffusivity of the ions beyond that of 
neutral molecules. Increased ion diffusivity results in a 
broader time profile and shifts the measured residence time ot 
the ion to times longer than those of the neutral molecule. 
Under the flow conditions of these experiments the CIMS 
instrument sampled about 18% of the flow r tube area. An ideal 
laminar flow with no radial or axial diffusion would result in a 
flow rate, averaged over a central cross section of the flow 


Here, K = k\ + kd, where k\ is the pseudo first-order rate 
constant for loss of H 2 S0 4 onto aerosol and kd is the 
background decay constant for diffusion and wall losses of 
H 2 S0 4 . The variable k n is the rate constant for reevaporation 
of H 2 S0 4 from the aerosol surface, and [H 2 S0 4 ] s is the aerosol 
surface concentration of H 2 S0 4 . 

Previous measurements by Marti et al. [1997] of the 
H 2 S0 4 vapor pressure above (NH 4 ) 2 S0 4 aerosols indicate 
that the evaporation rate of H 2 S0 4 from the aerosol is quite 
small, producing an H 2 S0 4 concentration of the order of 10 4 
molecules cm' 3 , below detection limit of these experiments. 
H 2 S0 4 coverage of the aerosol surface over the duration of the 
longest decay times was less than 1% of the total aerosol 
surface area. Because of the low H 2 S0 4 surface coverage and 
low H 2 S0 4 vapor pressure above (NH 4 ) 2 S0 4 aerosols, 
reevaporation of H 2 S0 4 from both (NH 4 ) 2 S0 4 and NaCI 
aerosols was considered negligible, and the second term in 
equation (2) was set equal to zero. 

A determination of k\ will depend on the mass transport 
between H 2 S0 4 and the particles. The particles generated in 
these experiments are between 20 and 400 nm in diameter and 
fall in the transition regime between gas kinetic and diffusion 
limits for mass transport. For uptake onto a monodisperse 
aerosol the first-order rate constant can be calculated in terms 
of a mass accommodation coefficient, a, with the Fuchs- 
Sutugin equation [Fuchs and Sutugin, 1970] 

<xN n cnD 1 

k { = £ (3) 

1 4(1 + afl) 


where 


(0.75 4-0,283^) 
Km Kn + 1) 


( 4 ) 


From the above, N p is the particle number density, c is the mean 
molecular speed, D is the particle diameter, and Kn is the 
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Knudsen number, which is a function of the molecular 
diffusivity, D g , and particle size: Kn - 6Dg/cD. Here the mass 
accommodation coefficient represents the number of H 2 S0 4 
molecules reversibly adsorbed onto the particle surface per gas 
kinetic collision. 

For a polydisperse aerosol distribution the observed k\ is 
the integral with respect to the particle diameter, D, of the size 
dependent pseudo first-order rate constant for H 2 S0 4 loss on 
the aerosol, i.e., 

k { =]k^D)dD (5) 

Numerically integrating over D, we obtain the following 
equation: 


acnDf 


dN 


* 1 = 1 - 


d\nD, 


40 + aflj) 


-AlnD, 


( 6 ) 


The aerosol number concentration for each particle diameter, 
( dN/d\nD)A\nD , was calculated as an average of all of the size 
distributions recorded at various injector positions for a 
given decay. Calculation of the mass accommodation 
coefficient involved iteratively changing a in the equation for 
k\ until the calculated value of k\ equaled the measured value. 

The diffusion coefficient, D g , of H2SO4 in N 2 was chosen to 
be 0.1 ± 0.02 cm 2 s' 1 in accordance with prior calculations 
[ Roedel , 1979; Marti et al. , 1997] and adjusted for an ambient 
pressure of 615 torr. In a humid airflow a fraction of the H 2 S0 4 
molecules are hydrated and have a lower effective diffusivity. 
However, at the low relative humidities of these experiments, 
most of the H 2 S0 4 molecules were unhydrated or 
monohydrated. The small decrease in diffusivity of the 
monohydrate relative to H 2 S04 is expected to be well within 
the uncertainty in the value of D g [Marti et al ., 1997]. 


Results 

The first-order decay of H 2 S0 4 onto an aerosol surface was 
measured for submicron polydisperse particles composed of 
(NH 4 ) 2 S0 4 , NaCl, and NaCl coated with stearic acid. Figure 2 
shows the H 2 S0 4 signal versus injector position for different 
concentrations of NaCl aerosol. The measured first-order rate 
constants and calculated values of the mass accommodation 
coefficient along with the uncertainties are given in Table 1 for 
H 2 S0 4 loss onto each of the substrates. The mass 
accommodation coefficients for H 2 S0 4 loss onto (NH 4 ) 2 S0 4 
and NaCl aerosol were found to be 0.73 + 0.21 and 0.79 + 
0.23, respectively. The measurement uncertainties are given as 
2a, twice the standard error for each set of measurements. The 
measured mass accommodation coefficients for H 2 S0 4 uptake 
onto NaCl particles coated with stearic acid had an average 
value of 0.25, indicating a significantly lower loss rate for 
HtS0 4 onto an aerosol with a hydrocarbon film. 

The mass accommodation coefficients measured for H 2 S0 4 
loss onto NaCl and (NH 4 ) 2 S0 4 aerosol in this study are 
similar and fall within the uncertainty range of either measured 
value. The relative humidity of the flow gas in these 
experiments was well below the deliquescence point of both 
salts, so that the uptake likely represents a physical 
adsorption of H 2 S0 4 onto a partially hydrated crystalline salt 
surface [Ewing and Peters , 1997]. The similarity in the two 
surfaces with respect to H 2 S0 4 adsorption is evident from the 
small difference in a for the salts. 

The value of a found for a NaCl aerosol coated with 
stearic acid was approximately a factor of 3 smaller than that 
measured for the bare aerosol. Stearic acid is a long-chain 
carboxylic acid (CH 3 (CH 2 )i 6 COOH) that is insoluble in 
water and consequently is expected to form a film on the 
aerosol surface. Figure 3 shows the aerosol size distribution 
for bare NaCl particles and NaCl particles with thick and thin 
coatings of stearic acid. The aerosol with the higher coverage 
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Figure 2. Relative H->S0 4 signal versus injector position under conditions of (a) no aerosol present, (b) an 
aerosol surface area of 3.0 x 1 O’ 6 cm 2 ern^, and (c) an aerosol surface area of 8.8 x 1 0 6 cm- cm L 
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Table 1. Summary of Measured First-Order Rate Constants (*i), Diffusional and Wall Loss 
Constant (kj), and Mass Accommodation Coefficients (a) 


RH, % 

M ' 1 

*d' s ' 1 

a 

Uncertainty 



NaCI 



1 6 

0.378 

0.120 

0.62 

0.20 

0 3 

0.139 

0.093 

0.85 

0.30 

1 7 

0.136 

0.108 

0.78 

0.25 

1 4 

0.169 

0.098 

0.91 

0.34 

1 1 0 

0.391 

0.048 

0.91 

0.26 

5.2 

0.341 

0.066 

0.70 

0.21 



(NH 4 ) 2 S0 4 



1.9 

0.383 

0.109 

0.75 

0.22 

5.2 

0.327 

0.066 

0.80 

0.26 

0.5 

0.250 

0.122 

0.64 

0.37 



Stearic Acid 



5.4 

0.241 

0.074 

0.31 

0.08 

4.4 

0.086 

0.089 

0.19 

0.06 


RH, relative humidity. 


of stearic acid resulted in a faster loss rate of Ff 2 S0 4 and thus a 
larger value for ot;(X = 0.31 and 0.19 for high and low stearic 
acid coverage, respectively. This difference between high and 
low coverage could be due to a change in the orientation of 
the acid on the aerosol surface, leading to a concomitant 
change in the number of polar binding sites available for 
H 2 S0 4 . Alternatively, the aerosol surface with the higher 
stearic acid coverage may have been highly irregular, resulting 
in a higher aerosol surface area. In the high-coverage 
experiments the appearance of a shoulder at large particle 
diameters in the particle size distribution may indicate the 
presence of stearic acid aerosols in addition to the coated salt 
aerosols. The shape of stearic acid aerosols may be even more 
irregular than that of the coated aerosols and can lead to an 
underestimate of the total aerosol surface area. Certainly, more 
knowledge of the aerosol surface structure and a far more 


extensive set of uptake measurements are needed to understand 
H 2 S0 4 loss rate onto aerosols containing organics. 

Discussion of Error 

Error in the calculated value of a will arise from 

uncertainty in the diffusion coefficient, point to point 
variation in the aerosol surface area during a single decay, 
uncertainty in the measured aerosol surface area, and changes 
in the H 2 S0 4 vapor pressure with instabilities in the flow gas 
relative humidity. Changes in the value of D g of ±20% lead to 
only a 3% change in a. An uncertainty in a near 5% could be 
accounted for from the standard deviation in the averaged 
aerosol size distributions within a single experimental decay 
measurement. In contrast, the FI 2 S 04 vapor pressure above the 
solution source was very sensitive to changes in the relative 


30000 ~r 



Figure 3. NaCI aerosol size distributions with no stearic acid (triangles), low stearic acid (diamonds), and 
high stearic acid coverage (circles). 
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humidity, which added to the uncertainty in the observed 
first-order rate constant and mass accommodation coefficient. 
The vapor pressure of H 2 S0 4 above a solution declines 
exponentially with increased relative humidity [Marti et ai, 
1997; Ayers et ai , 1980], In the dry flow gases of these 
experiments, a 0.1% change in the relative humidity (RH) in 
the flow gas from 1 to 1.1% could produce a change in the 
H2SO4 vapor pressure of the order of 3 x 10 7 molecules cm' 3 , 
i.e., close to 15% of the concentration range of the signal 
decay. The H 2 S0 4 signal was far less sensitive to changes in 
the water vapor pressure at a relative humidity above 5%. 
Error arising from small variations in the flow gas relative 
humidity and aerosol size distribution is evident in the 
calculated fit of the signal to an exponential decay. As can be 
seen from Table 1, the uncertainties in the calculated values of 
a are highest in the decays with the lowest relative humidity. 

As a result of laminar flow and diffusion the H2SO4 
concentration will develop a parabolic profile across the flow 
tube. Because of this radial gradient in H 2 S0 4 the observed 
decay constant is a nonlinear function of the first-order wall 
loss rate and the pseudo first -order aerosol loss rate. To 
quantify the error introduced by our simplified analysis 
(which assumes that the observed rate constant is a linear sum 
of k\ and kj), we determined k\ using the method of Brown 
[1978], modified for species detection in a central cross 
section of the laminar flow. With this method we calculated a 
cross-sectional sampling area based on the total flow through 
the sample inlet and the laminar flow velocity profile in the 
flow tube. Using this method, we found calculated values of 
first-order rate constant, k\, within 2% of those values 
calculated by using the measured flow time from equation (1). 
The standard errors of the average value of a in both methods 
were the same. Here the radial gradients in the flow time and 
H 2 S0 4 concentration across the sampled region were small 
enough that the obsen/ed rate constant, K, could be 
considered a linear sum of k\ and kd. 

Error in the value of a will arise also from uncertainty in 
the SEMS measurements of the aerosol size distribution and 
calculated particle surface area. Most of this uncertainty is 
associated with accuracy in the pressure, flow rates, dilution 
of the sample gas, transmission of the particles through the 
DMA, the particle charging probabilities, and the inversion 
algorithm used to characterize the particle size distribution. 
We estimate the sum of the uncertainties in these parameters to 
be approximately 25% of the calculated particle surface area. 
The reported uncertainties in the average values of a for loss of 
H 2 S0 4 onto NaCl and (NH 4 ) 2 S0 4 particles is a sum of the 
25% systematic error associated with the uncertainty in the 
aerosol surface area, a 3% error due to uncertainty in the 
HtS0 4 diffusion coefficient, a 3% error in the flow' time 
measurement, and twice the standard error (2a). 


Atmospheric Implications 

Under the ambient conditions of the lower troposphere the 
composition and phase of submicron particles depend on the 
local sources of condensable gases and the relative humidity. 
The fine-mode particles in a remote marine environment are 
primarily aqueous solutions of sulfate, methyl sulfonate, and 
ammonium ions [Heubert et ai, 1996; Savoie et ai , 1993]. In 
addition to sulfate and ammonium ions, aerosols in 
continental air masses may contain nitrate ions and as much as 


20-50% by mass organic molecules [Novakov and Penner , 
1993; Saxena and Hildemann , 1996]. Also, gas phase sulfuric 
acid will be hydrated with one or more H 2 0 molecules almost 
entirely throughout the lower troposphere [Kulmala et ai , 
1991; Marti et ai , 1997]. 

We can expect the uptake rate of H 2 S0 4 onto particles to 
vary with the level of H 2 S0 4 hydration and the particle 
composition. Although these experiments were conducted 
well below the deliquescence point of NaCl (75% RH) and 
(NH 4 ) 2 S0 4 (80% RH), the aerosol surfaces likely had some 
H ? 0 coverage. Studies of H 2 0 adsorption isotherms on NaCl 
crystallites by Barraclough and Hall [1974] indicate that a 
NaCl crystallite surface has <0.1 to 0.5 monolayer of H 2 0 
coverage over the corresponding 0.3-11% relative humidity 
range of the present experiments. Water adsorbed to the 
aerosol surface may provide a polar binding site and thereby 
enhance H 2 S0 4 uptake. Also dependent on relative humidity 
is the percent hydration of H 2 S0 4 . At a relative humidity 
below 1%, H 2 S0 4 is almost entirely unhydrated. However, at 
a 10% relative humidity, approximately 50% of the gas phase 
H 2 S0 4 will be in a monohydrate form [Marti et ai y 1997]. No 
appreciable trend in the mass accommodation coefficient with 
relative humidity was observed over the limited range of 
0.3-11% RH. The small difference in the mass accommodation 
coefficients of H 2 S0 4 for (NH 4 ) 2 S0 4 and NaCl particle 
surfaces of the present study indicates either a similarity in the 
degree of H 2 0 coverage on these two salt particles or the fact 
that the uptake of H 2 S0 4 onto the aerosols was relatively 
indifferent to H 2 0 on the salt surface. In the case of H 2 S0 4 
uptake onto a neutral salt the binding energy of H 2 S0 4 
appears strong enough that the acid readily adsorbed to the 
partially hydrated salt surface, and further enhancement of the 
uptake rate by additional water coverage may have little effect 
on an already high mass accommodation coefficient. 

The uptake of H 2 S0 4 onto wet, ambient aerosol droplets 
involves the adsorption and subsequent solvation of H 2 S0 4 
into the aqueous droplet. Here the energetics of solvation and 
hence reevaporation rate of H 2 S0 4 depends on the Henrys law 
constant, which varies with the particle composition, pH, and 
temperature. The vapor pressure of H 2 S0 4 above a mixed 
sulfate aerosol (i.e., a sulfuric acid/ammonium sulfate solution) 
was found to increase with the percent H 2 S0 4 composition of 
the aerosol [Marti et ai , 1997]. Thus the net uptake of H 2 S0 4 
is expected to decrease with particle acidity. 

As is evident in the observed decrease in a for an aerosol 
coated with stearic acid, the initial adsorption of H 2 S0 4 
varies with the aerosol surface composition. As an oxidizing 
agent, H 2 S0 4 may react chemically with organics present in 
the aerosol, leading to increased adsorption of H 2 S0 4 onto 
aerosols. Our initial attempt at coating the aerosols with 
stearic acid involved solvating the stearic acid in ethanol or 
propanol. However, we found that these alcohols readily 
reacted with H 2 S0 4 to form organic-sulfate clusters, which 
completely depleted the HS0 4 ‘ signal. This serendipitous 
discovery may indicate that H 2 S0 4 reacts with atmospheric 
organics, either in the gas phase or in aerosol solutions. 

The measured mass accommodation coefficients for NaCl 
and (NH 4 ) 2 S0 4 aerosol may coincide well with the loss rate of 
HoS 0 4 onto neutral salt aerosols found in a remote marine 
environment. However, evidence from the reduced loss rate oi 
H,S0 4 onto an aerosol coated with an organic film and a 
previous experiment that found the H 2 SG 4 vapor pressure 
above an aerosol to increase with aerosol acidity [ Marti et ai , 
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1997] suggests that H 2 S0 4 uptake onto aerosols varies 
broadly with the particle composition. Field measurements at a 
remote continental site that measured steady state H 2 S0 4 
production and loss rates found good agreement with their 
data, using a value of a of 0.5 and 1.0 [Eisele and Tanner , 
1993; Weber et ai , 1997]. In these two Field studies the 
aerosol likely had a mixed aqueous composition of S0 4 2 -, 
NO 3 ', NH 4 + , inorganic metal and organic ions. 

Conclusion 

In this study we measured directly the H 2 S0 4 loss rate 
onto submicron aerosols for three different aerosol substrates. 
The mass accommodation coefficients were calculated from the 
first-order rate constants and were found to be 0.73 + 0.21 and 
0.79 + 0.23 for H 2 S0 4 loss onto (NH) 4 S0 4 and NaCl aerosols, 
respectively. The loss rate of H 2 S0 4 onto a NaCl aerosol 
coated with stearic acid was lower than that of the salt aerosol 
with measured values of the mass accommodation coefficient of 
0.31 and 0.19 for high and low stearic acid coverage, 
respectively. 

The values of a are considerably higher than previous 
measurements of H 2 S0 4 uptake onto aqueous H 2 S0 4 droplets 
[Van Dingenen and Raes , 1991] but compare well with the 
value of a = 0.5 to 1.0 deduced from the results of a field study 
of steady state production and loss rates of H 2 S0 4 at a remote 
continental site [ Eisele and Tanner, 1993; Weber et al 
1997], The higher values of a will result in higher calculated 
aerosol growth rates. In steady state model calculations of 
aerosol formation and growth, a higher aerosol growth rate 
would mean lower estimates of the gas phase H 2 S0 4 
concentration and consequently slower aerosol nucleation 
rates. 

The uptake of H 2 S0 4 was found to depend strongly on the 
aerosol composition, with a lower value of a reported for an 
aerosol with an organic coating. A similar reduction in a has 
been observed for uptake of NH 3 onto aerosol with an 
organic coating [Daumer et al 1992]. The variability in a 
with aerosol composition suggests that model evaluations of 
aerosol nucleation and growth need to consider the regional 
aerosol composition. In their model study, Kerminen and 
Wexler [1995] proposed using a value of a that decreases 
with increased aerosol size, because larger particles may have 
a higher organic composition and therefore a slower H 2 S0 4 
uptake. However, the chemical reactivity of H 2 S0 4 with 
organics may affect their uptake onto aerosols such that a may 
vary with the specific organic compounds present in the 
aerosol. 

Further studies of H 2 S0 4 uptake onto aerosols indicative 
of the lower troposphere are needed to improve our 
understanding of aerosol growth and H 2 S0 4 loss rates. 
Studies of the uptake rate of H 2 S0 4 onto an aerosol surface of 
varying organic composition as well as aerosol pH, 
temperature, and relative humidity would enhance the ability 
of model calculations to predict aerosol formation and growth 
rates over a range of diverse atmospheres. 
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